Casimir-Lifshitz force out of thermal equilibrium and asymptotic non-additivity 
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We investigate the force acting between two parallel plates held at different temperatures. The 
force reproduces, as limiting cases, the well known Casimir-Lifshitz surface-surface force at thermal 
equilibrium and the surface-atom force out of thermal equilibrium recently derived by M. Antezza 
et at, Phys. Rev. Lett. 95, 113202 (2005). The asymptotic behavior of the force at large distances 
is explicitly discussed. In particular when one of the two bodies is a rarefied gas the force is not 
additive, being proportional to the square root of the density. Nontrivial cross-over regions at large 
distances are also identified. 
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The study of the thermal fluctuations of the electro- 
magnetic field and of their effects on the force acting on 
surfaces and atoms is a longstanding subject of theoret- 
ical research starting from the seminal Lifshitz paper [l| 
(see also Q). The elusive nature of the thermal compo- 
nent of the force follows from the fact that thermal ef- 
fects becomes visible only at large distances, of the order 
of the photon thermal wave length At = hc/ksT, where 
they prevail on the Casimir force originating from the 
T = quantum fluctuations of the field. At room tem- 
perature the thermal wave length corresponds to about 
7 microns, a distance at which both the Casimir and the 
thermal forces are very weak and difficult to reveal ex- 
perimentally. The existence of thermal effects has been 
experimentally demonstrated only recently by the JILA 
experiment [3|, by measuring the frequency shift of the 
center of mass motion of an ultracold atomic cloud lo- 
cated at a distance of a few microns from a dielectric 
substrate 0, [1]. 

Thermal fluctuations determine the asymptotic, large 
distance behavior of the electromagnetic pressure which 
takes the Lifshitz form |2] 
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in the case of two parallel surfaces separated by a dis- 
tance I. Here eig = ei(0) and £20 = £2(0) are the static 
dielectric constants of the two materials and T is the 
temperature of the system. When the system is not in 
thermal equilibrium the pressure is expected to exhibit 
a different behavior. In particular in recent papers the 
Trento team 0, § has shown that the surface-atom 
force out of thermal equilibrium exhibits a new asymp- 
totic behavior at large distances. With respect to that at 
equilibrium, the new force exhibits a slower dependence 
on the distance and a stronger dependence on the temper- 
ature, making its experimental detection more accessible 
as demonstrated in the experiment of 

The purpose of the present work is to investigate the 



behavior of the force out of thermal equilibrium in the 
case of two parallel surfaces. The general goal is to bet- 
ter understand the role of thermal fluctuations which is 
not fully exploited at thermal equilibrium, being masked 
by peculiar cancellation effects between propagating and 
evanescent waves [1, 01- In particular we address the 
following questions: a) how is the Lifshitz law - and its 
asymptotic limit ([T]) - modified if the temperatures of the 
two bodies are different? b) Can one recover the results 
of for the surface-atom force when one body is made of 
a very dilute material corresponding to a gaseous phase 
with (e - 1) ^ ? 

Let us consider two parallel dielectric half spaces locally 
at thermal equilibrium with different temperatures and 
separated by a distance I. In our configuration the left- 
side (right-side) body has a complex dielectric function 
ei{uj) (£2(w)) and is held at temperature Ti (T2), the 
whole system being in a stationary configuration. We 
assume that each body fills an infinite half-space (see, 
however, discussion after Ea. (|17p ). In practice this means 
that the bodies are thick compared to the penetration 
depth of the thermal radiation. In such conditions the 
presence of the remote surfaces of the bodies results only 
in a ^-independent contribution to pressure, which will 
not be considered in this Letter. 

Let us assume that the separation between the 
bodies is in the z— direction. Then the electromag- 
netic pressure between them is given by the aver- 
age P"<'q(Ti,r2,0 = (T,.(r,t)), where T,,{v,t) = 
{El - El- Bl- Bl- Bl) /Stt is the zz compo- 
nent of the Maxwell stress tensor Q in the vacuum gap. 
In this work we focus on the thermal component Pth of 
the force, defined by P(Ti, T2, /) = Po(0 + ^'th(Ti, T2, /), 
where Pq is the T — Q quantum pressure originating from 
the vacuum fluctuations of the field 0]- 

The electromagnetic field in the vacuum gap physi- 
cally originates 9] from the fluctuating polarization fields 
P[ti;;r] in the two half-spaces whose correlations inside 
each body are described by the fluctuation-dissipation 
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FIG. 1: Relative contribution of the thermal component with 
respect to the zero-temperature component of the pressure 
between two different materials: Fused Silica (Si02, body 1) 
and Silicon (Si, body 2). 



spaces add incoherently. Assumption ([2]) (local source 
hypothesis) was first used in and represents the start- 
ing point of our analysis allowing for an explicit calcula- 
tion of the electromagnetic field also if the system is not 
in global thermal equilibrium The electric 

field at the point r in the gap can be in fact expressed in 
terms of the source polarization field via the convolution 
Ejw; r] = / G [w; r, r'] • P [oj; r'] dr' of the Green tensor 
[l4| , where the integration is performed over the volume 
containing the sources at r'. At the same time the mag- 
netic field is easily evaluated using the Maxwell equation 
B[uj; r] = -iV A E[uj; r]/k, where k = uj/c. Then Eq. ^ 
allows us to write the thermal pressure acting between 
the bodies as the sum of two terms: 

priTi,T2,i) = nrm, 0, o + nr(o, t^, o, (s) 

each of them corresponding to a configuration where only 
one of the two materials is at non-zero temperature. It 
is convenient to write 
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where 2('^) is the imaginary part of the dielectric func- 
tion of the two materials. Due to the presence of the 
i5(r — r') factor these correlations are local so that the 
effects of the fluctuations originating from the two half- 
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where P^^{T,l) is the Lifshitz pressure at equilibrium 
2]. If we write the electromagnetic wave vector in the 
vacuum gap as k = (Qx, Qy, ^z), whose longitudinal part 

has modulus Q = ^ Q"^ + Qy, and Qz = \/k'^ — , it is 

possible to express 1;he quantity APth as the sum APth = 



APPW 



I^Pf^ , with 



1 

,hu)/kBT _ I 



dQ Q q, 



E 



1 - Ir 



1 '2 I 



APEW(T,0 



^ rd.. ^ r^nn, -2nn... Im(rnRe(r,^) - Im(r^)Re«) 



(6) 



(7) 



where we have separated the effect of propagating waves 
(PW) from that of evanescent waves (EW) correspond- 
ing, respectively, to real and imaginary values of Qz. We 
have also subtracted the /-independent contributions so 
that both Eqs.([6]) and ([7]) vanish as I oo. In Eqs. 
(ini) and ([7|) the dielectric properties of the two materials 
enters through the reflection Fresnel coefficients for the 
vacuum-dielectric interfaces 

(m) (m) 
s _ Qz ^ Qz p _ 9z£m ~ Qz 

^™ ~ ', M ' ~ . (rn)' ^ ' 

qz + Qz qz£m + ql 

where s and p correspond to the transverse electric and 
magnetic polarizations and qi'"'* = a/ Em k"^ ~ is the z- 
th component of the wave vector in the material m = 1,2. 



The effect of multiple reflections between bodies is ac- 
counted for by the denominator = 1 — r'^r'^e^^'^^K It 
is evident that the pressure P^^'^{Q,T,l) of Eq. (O can 
be obtained from Eq. ^ by replacing <-> rj. It is 
also immediate to see that Eq. (jSj at thermal equilib- 
rium Ti = T2 = T provides the well known equilibrium 
pressure P^^[T, I). For identical bodies = r!^, yielding 
APth — 0, the pressure ^ is given by the expression 

p,rm, = p:^iTi,i)/2 + p:^{T2,i)/2. (9) 

Equation previously obtained in [l^, is not however 
valid if the two materials are different. This is in dis- 
agreement with the results of [1^1, where @ was found 
to be valid in general [l3|. In particular Eqs. ^ and 
([5]) show that interchanging the temperatures of the two 
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materials implies a change in the force if the two materi- 
als are different. This clearly emerges from FiglU where 
the results of a full calculation of the pressure between 
Fused Silica (body 1) and Silicon (body 2) are presented 
[ri |. It is interesting to note that the relative values of 
the pressure for the two non-equilibrium configurations 
(dashed and solid lines) strongly depends on the temper- 
atures of the two bodies, on the separation and on the 
positions of the resonances of the two dielectric functions 
£i^2{^)- One can also see that both values are smaller 
than the one at thermal equilibrium (dotted line). This 



it is not however always the case. In fact if one of the 
two bodies is rarefied the non equilibrium pressure can 
become larger than at equilibrium (see and discussion 
below) . 

In the following we will focus on the behavior of the 
force at large distances. For this study we will consider 
only dielectric bodies (where the limit £1,2 (w — > 0) is 
finite), while the case of metallic bodies will be presented 
elsewhere. The expansion of Eq.(|3]) results in the large 
distance behaviors 
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holding for I ^ At (if one of the bodies is rarefied the 
condition becomes more stringent, see Eq. lfTil) below). 
Here C(3) is the Riemann function, r^(i) are the Fresnel 
reflection coefficients calculated from Eq.® setting the 
static approximation e„i = £.mo, and = (l -l-t^). 
The PW and EW contributions of Eqs. ^ and HID 
are of the same order but have opposite signs, the PW 
term being attractive and the EW repulsive. This feature 
is not however general and at shorter distance the signs of 
the PW and EW terms can change as discussed below. At 
distances of the order of the thermal wavelength Eqs.([TU]) 
and pT|) , as well as the Lifshitz result H]) , provide only a 
crude estimate of the pressure. For example in the case 
of the Fused SiHca (£10 ~ 3.8) - Silicon (£20 ~ 11.7) con- 
figuration at I = 5/im, the asymptotic laws overestimate 
the full calculation shown in Fig.(IT|) by a factor 1.5, 2 
and 3 for the dashed, dotted and solid lines, respectively. 

Since in this work we are also interested in recovering 
the surface-atom force which is relevant for the recent 
experiments with ultracold gases 3 , it is useful to study 
the case in which the body 2 is cold and very rarefied, 
corresponding to small values of (£2 — 1) = Aima. Here 
n is the density of the material 2 and a is the dipole 
polarizability of the atoms. The expansion of pressure 
in Eq. Q should be performed through two limiting 
procedures: the large distance I — s- 00 and the diluteness 
(£2 — 1) conditions. It is crucial to identify the proper 
order of the two limits. 

One relevant asymptotic behavior is obtained by first 
taking the limit of large I at fixed £2 (this yields, by the 
way, expressions (|10l) and (|lip ). and then carrying out 
the limit of rarefied body. One finds the expression 
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for the total pressure, where C — Cpw + Cew 
3.83 • 10~^, is a numerical factor with Cpw — C(3)/87r 
4.78 • 10-2 and Cew ~ -0.96 • lO'^. 



The peculiar -\/£20 — 1 dependence of p2)) means that 
the pressure acting on the atoms of the substrate 2 is 
not additive. Additivity would in fact require a linear 
dependence on (£20 — 1), and hence on the gas density n, 
as happens for the Lifshitz pressure ([1]) as (£20 — 1) 0: 
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The non additivity of the pressure follows from the 
fact that for large I the main contribution in the force is 
produced by the grazing waves incident on the interface 
of the material 2 from the vacuum gap at small angles 
corresponding to small values of g^/fc < ^£20 — 1- Hence 
the reflection coefficients from the body 2 is not small 
even at small £20 — 1 and the body cannot be considered 
dilute from the electrodynamic point of view. This is a 
peculiarity of the non-equilibrium situation. In fact at 
equilibrium this anomalous contribution is canceled by 
the waves impinging the interface from the interior of the 
dielectric 2, close to the angle of total reflection. Notice 
that result IT^ is valid at the condition 



I > At/V£20 - 1, 



(14) 



which becomes stronger and stronger as (£20 — 1) — * 0. 

The second limiting procedure is obtained by first tak- 
ing the expansion of ([4]) for small values of (£2 — 1), 
and then carrying out the limit of large distances. In 
this case the relevant wavevectors satisfy the condition 
^ \/£20 — 1, and the PW component produces a 
contribution identical to the EW one, yielding 
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in units of the asymptotic behavior (|T2| 



FIG. 2: Dimensionless function f{x) (see Ea. (fT7|l ) 
dashed line gives the asymptotic limit at small x. 
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15)) holds in the distance range complementary 



At < ^ < Xt/V£20 - 1- 



(16) 



In deriving Eq. ()15p we also replaced £1,2 (w) with their 
static values Emo- This is justified if fc^T is much smaller 
than the lowest resonances of both the body 1 and the 
atoms of the dilute body 2. It is worth noting that the 
interval (|16p practically disappears for dense dielectrics. 

It is also worth noticing that, due to the diluteness 
condition (£20 — 1) <C 1 and as a consequence of Eqs.([4|) 
- dS]), in both regions p4|) and (fTH]) the term APth gives 
the leading contribution to the pressure and consequently 
the dependent interaction between the two bodies will 
be attractive if T\ > T2 and repulsive in the opposite 
case. 

The transition between the two regimes ((T^) and p5|) 
can be investigated performing the diluteness limit (£20 — 
1) ^ in Eq.ijl]), for a fixed value of the dimensionless 
variable x — l\/e2o ~ 1/At- The results are reported in 
Figl2l where the thermal pressure P^^'^iT, 0, 1) is plotted 



p;r(T,o,o 
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Here f{x) is a dimensionless function of the variable x. 
When a; 00 (regime ([12])) one has f{x) 1, while 
when X — > (regime (|15l) ) one find f{x) x/2AC w 
1.09a:. 

In order to recover the asymptotic result of |6| for the 
surface-atom force out of thermal equilibrium it is crucial 
to follow the second limiting procedure, leading to result 
P^ . In this case, however, the PW term must be omitted 
since the atomic gas occupies a finite region of space and 
does not absorb the thermal radiation. Using the formal- 
ism of the present work this corresponds to treating the 
body 2 as a slab of rarefied gas of thickness L for which 
one should also take into account the force acting on its 
remote surface. In the absence of absorption [l9| it is 
possible to show that, including refraction at the remote 
surface, the PW pressure becomes vanishingly small, of 
order (£20 ~ 1)^ with respect to Eq. ([15]). In this case 
the EW contribution, which is 1/2 of (fT5|) . provides the 
total pressure actingon the gas and is exactly equivalent 
to equation (12) of @ for the surface-atom force. Notice 
that in the derivation of Q the leading role of the EW 
term was stressed from the very beginning. 

In conclusion we have generalized the Casimir-Lifshitz 
surface-surface force to configurations out of thermal 
equilibrium and calculated the corresponding asymptotic 
behavior. When one of the two bodies is a rarefied gas 
a cross-over region emerges where the pressure changes 
from a T^/P behavior, characterizing the surface-atom 
interaction, to a region at very large distances where the 
pressure behaves like T/l^ and is no longer additive. 
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